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We study the electric dipole absorption in small metal particles and thin films in a longitudnal 
electric field. In diffusive approximation, we give both the phenomenological and microscopic deriva- 
tions with the account for Thomas-Fermi screening and Drude relaxation. 

I. INTRODUCTION 

Classical absorption in a variable electromagnetic field is a well studied and documented field of resear ch §. Here, 
we will be primarily concerned with the absorption of the longitudinal component of a "quasi-static" electric field 
EQe lwt . This problem arises, for instance, in study of absorption of electromagnetic radiation by particles whose size 
is smaller than the wave-length A [0, |§j or for absorption in a film sandwiched between the capacitor plates. We 
point out the role of the diffusion, in addition to ohmic conduction, in the dynamics of the charge transport. 

Absorption is determined by the imaginary (" dynamical" ) part of the electric moment Pe luJt |l|] 

Q = ^Re{iuVE } = -~uE Im{V} (1) 
For metals, the standard approach is to characterize the particles by the Drude dielectric constant, 

e D = l-^ (2) 



where 



°d = -7—. — (3 

1 + VjJT 



is the Drude conductivity and oo and r are the Boltzmann conductivity and the scattering time respectively. 

At small frequencies, er> is dominated by the imaginary part. Consequently, that the quasi-static polarization 
P s t ~ Eo/Att, attributed to the surface charge, screens the external field and produces the electric moment 

V st = P st V ~ M 
An 

where V is the volume. The dynamical polarization P u ~ —Eq/Attcd produces the electric moment 

iuE V 



(Air) 2 a 



which, upon substitution in eq. (Q), yields the standard Rayleigh-Drude absorption 

(LuEpfv 
2 (4tt) 2 ctq 

Another way to interpret this result is to notice that the current density in the bulk is given by 

. n iuE 
~ uvPst : 

47T 
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This current is responsible for the surface charge that screens the external field. It couples to the dynamical field in 
the bulk 



resulting in absorption 

The surface charge is, of course, only a mathematical construction. In reality, the field is screened over the screening 
length A -1 , where A is the Thomas-Fermi wave- vector 

A 2 a 2 47TC70 

A = 47re 2 — = — — - (4) 
dp D w 

dn/dp, being the thermodynamic density of states and D is the diffusion constant. For a good metal A -1 is very small 
- typically an Angstrom - and the sample size L is such that LA 3> 1 . Furthermore, tt> 1 where I — \J Dt is the 
electron mean-free-path. In semiconductors, on the other hand, various relationships between L, I and A are possible. 
The Drude dielectric constant obviously cannot describe the charge distribution of the screening layer. This paper 
presents a natural extension of the classical argument based on o\d alone to incorporate both static and dynamical 
aspects of the finite thickness of the screening layer. 

II. PHENOMENOLOGICAL DESCRIPTION OF ABSORPTION 

A. Maxwell equations and boundary conditions 

Consider a sample in the quasi-static field E z = Eoe luJt . Using the Coulomb gauge for the longitudinal field, we 
solve together the Maxwell equation and the continuity condition 

47175(1-; = V-E(r;w) (5) 
-iwp(r;u) = V- j (r;w) (6) 



Eliminating p, we find 



V- (j(r;w)+^E(r; W ) ) = 



whereof the current density can be written as 

j(r; W ) = ^(D(r;a,)-E(r;u;)) (7) 

in terms of the divergenceless vector D 

V-D(r;w) = (8) 

The boundary conditions complementing eq. (|7]) can be obtained as follows. First, there is no current through the 
sample boundary d 

j n (r;w)\ d = (9) 

and, consequently, 

D n {r-u))\ d = E n (r;u)\ a (10) 
Second, since p is bounded and there is no surface charge, the electric field must be continuous at the boundary 

E(r; W )| 9 = E°"*(r;a,)| a (11) 
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where E OM * is the held outside the sample. 

D has a meaning of electric displacement held for a " dielectric" medium in which all the current is assigned to the 
polarization P, namely, 

j(r;w)=iwP(r;w) (12) 

so that 

D(r;w) = E(r;w) + 4nP(r;uj) (13) 

and 

V.P(r;w)=-/j(r;w) (14) 

Eq. (0) can be used also to dehne T) out — E ou * outside the body. From eqs. ( pd] ) and (|l(]) it then follows that the 
normal component of D would be continuous at the boundary. On the other hand, the tangential component of D 
may be discontinuous due to the tangential current at the boundary. 



B. Calculation of the absorption 

Local absorption is given by 

Q (r; uj) = ~ Re {j (r; u) -E* (r; u)} = - ~ Im {ujP (r; w) -E* (r; u)} (15) 
Using eq. (R), the latter can be also written as 

Q(r;o;) = -Im{|^D(r;u;)-E*(r;a;)} (16) 
It must be emphasized that in evaluation of the total absorption Q 

Q = I Q{v-uj)dv (17) 



the actual field in the interior of the sample can be replaced by the applied field (i.e. in the absence of the sample). 
For instance, using eqs. (||) and (10), E* (v\uj) = — V<p* (r;w) and the Gauss theorem, we hnd the total absorption as 

W Im{ [ p(r;u,)(t>*(r;u)dr\ (18) 



(19) 
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It is shown in Appendix A that the latter can be also written 



Q= |lm|y p(r;cu)<f>o (r) drj = ~ Im {e -^} 



where 

0o (r) = -E -r (20) 

is the potential of the applied held and 

~ft = J p(r;u;)rdr (21) 

is the electric moment of the system (of course, the second eq. ( |l9| ) also follows from a general thermodynamic 
argument §). The physical reason why either <f> or </>o can be used in the dissipation integral is that their difference is 
due to the screening charge itself and is the (non-dissipative) energy of interaction of the screening charge with itself. 
By means of the Gauss's theorem and eqs. ([l4|), (||) and (|l2|), we also find 

Eo-7? = J - [V- (P (r; w) (E -r)) - E -P (r; w)] dr = E - J P (r; w) dr 

that is 

"P= J P(r;w)dr (22) 
which is consistent with the fact that E can be replaced by Eq in the evaluation of Q. 
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C. Constitutive equation and solution for the field 



A constitutive equation is needed to complete the description of the medium. For this we use the "generalized 
Einstein transport equation" 

j(r;w) =a D E(r;u)-D D Vp(r;u;) (23) 

where 

Dd = (24) 

1 + ILOT 

Eq. ( p3[ ) will be further discussed in Section III. Combining eqs. ^) and (|2^), we find, with the help of eq. (pi) the 
following equation for the electric field 

V 2 E(r; W )-A 2 E(r;c) = --^D(r;c) (25) 



l2 = ^B =K 2f 1 + _^_\ (26) 



(27) 



Here A is the dynamical screening length given by 

iujCj 

Solution of eq. (|2fj) is found as 

E (r; W ) = ^i^ + E' l (r;c) 

(D 

where E h is the solution of the homogeneous equation 

V 2 E' 1 (r; - A 2 E' 1 (r; lj) = (28) 

such that the boundary condition ([ll]) is satisfied. Note that E' 1 represents the deviation from the standard Rayleigh- 
Drude description. 

The procedure for finding the electric field can be, in general, described as follows. First, one determines the 
displacement vector D = — V</> d , 

VV = (29) 
V n <f> d \ 9 = -E™ t \ d (30) 

as per eqs. (|^) and (|l0[), and the solution of the homogeneous equation 

Vy^rjaO-AVfaw) =0 (31) 

V n <f\B = -E^\e(^^) (32) 

as per eqs. (pSj), ( p7| ) and E h = —V<j) h . Subsequently, the total field in the interior is found as 

= <f> h + — (33) 
<P\d = -<p out \d (34) 

as per eqs. ( |27j ) and (|TT|). 

Clearly, the use of eq. ( p3[ ) assumes the applicability of the diffusion approximation. The conditions for the latter 
are ui < r _1 and, strictly speaking, £A <C 1. For £A 1, eq. ( p3| ) can still be used in the bulk. The situation within 
the screening layer at the boundary must be analyzed separately. 
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D. Examples: fields and absorption in a slab and sphere 



1. Solution for a slab 

Consider a narrow slab, L z = L <C L x , L y ~ ^/S. Transverse degrees of freedom can be neglected so the only 
variable is z and all fields are directed along z. The boundary condition is that the electric field inside and outside 
the slab match at the boundary. This is because the field of the slab's electric moment is zero, £ = 0. Indeed, this 
field has to satisfy the equation = (no charge) and fall off to zero at infinity (it is also just a field of two infinite 
plates of opposite charge). Thus D — E out = E and the boundary condition is given by 

E(0;u) = E{L;lo) = E (35) 

Eqs. (R) and d2q) become respectively 



and: 



dz 2 v ' ' v ' ' D D 

Solving eq. (pTl), subject to the boundary condition (pq), we find 



j(z;u) = -^(E(z;u>)-Eo) (36) 



4^E(z;lu)-A 2 E(z;lu) = -^E C!7) 




cosh ( (z - -§0 A 

E (z;uj) = | I + (e D 1) V L ~ ' \ — (38) 

cosh f A 

, , cosh ((z - 4) A) , s 

E st {z)=E V \ - T 2 / 1 (39) 



and 



cosh -j A 



L\ (e D -l) 



sinh 



((• 



L 




for the field and potential inside the slab, with the static components explicitly emphasized. 
Local absorption is given by 

Q (z) = ~ Re {j (z; to) E* (z; u)} = ^ljE Im {E (z; u)} (42) 

The total absorption of a slab is 

Q = s[ Q (z) dz = -j-wE S [ lm{E(z;ui)}dz (43) 
Jo 87r Jo 

= ~uE Im{V} (44) 



where the latter expression - in agreement with eq. (|19|) - is written in terms of the electric moment of the slab 

V = S [ zp(z;w)dz = ^-- -j-S [ E(z;u;)dz (45) 
J 4tt 4tt J 

and integration by parts has been used. Further integration of by parts yields 
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and combining this with eq. (pQ), we obtain the following expression for the electric moment of the slab: 



r _VE (e D -l) ( tanh(£A/2 ; 
4tt e D I LA/2 1 

^ = VEo L_ tanh(LA/2)X 



4tt V iA / 2 / 

where the static component is explicitly emphasized. Using eqs. ( [H] ) and (fl7|), we find the closed form solution for 
the total absorption in the slab 



Notice that no assumption about the magnitude of LA has been made. For LA 3> 1, expanding eq. (|49| ) at small 
frequencies, w <C D/L 2 , we find 

The expression in front of the parentheses is the standard Rayleigh-Drude result Q rd . The second term inside the 
parentheses is the screening layer volume effect. 



2. Solution for a sphere 

The field outside the sphere of radius a, E ont = — V0°"*, where Jl[ 



j nut 



7^- 



(r;u,) = -E -r+— (51) 



r 



and V || Eo is the sphere's electric moment. 4> out is proportional to the scalar Eo-r oc 003(6*) and so is the interior 
solution (j> d which also obeys the Laplace equation. Matching according to eq. (pG) gives 



^(r;c) = -E -r(l + -^) (52) 



Ea a 



Consequently, the field inside the sphere is 



E (l 



2V 



E(r;w) = * ^-+E h (r;w) (53) 



with the homogeneous solution found as E' 1 ~ ~V<fr h , where 

, flEn-r 



^(rjw) = A?=^i 1 (r&) (54) 
k (x) = ^J. (x) (55) 
and i\ (x) and 1 3 (x) is the spherical Bessel and Bessel function of imaginary argument respectively. The constants 



A and V are found from the boundary conditions (|ll|) (or, equivalently, eqs. (31) and (p3|)) and, using the notation 



coth I aA ) 

X = 1 - 3 + — ^ (56) 

UP 2 



aA 



G 



we obtain the following expression for A 
and the electric moment of the sphere: 



?>=^M^(^-1)] (57) 



4tt \e D + 2X 

3VE f 3 coth(aA) 

I aA ' ( a A) 2 



Vst = ^ ( 1 - + - (58) 



where the static component is explicitly emphasized. The latter coincides with the result of Ref. || obtained quantum- 
mechanically. 

Using eq. ([!]), we find the closed form solution for the absorption in a sphere 

Q = ±uE*VIm{( z X nv (tD-l))} (59) 



8tt u [\e D + 2X 

where, again, no assumption about the magnitude of a A has been made. For aA 3> 1, the small- frequency expansion, 

2(47r) 2 a V ^aA) \ ^Aj K ' 

If the oscillating field is due to electromagnetic radiation with wave-length longer than the sphere size, then dividing 
by the density of the incoming flux find for effective cross-section of the absorption 

4na c V 2aA / v ; 



III. RESPONSE-FUNCTION FORMALISM 



Using standard Fermi- liquid considerations Q in real space or a field-theoretical argument Q , it can be shown that 
the linear response of the current to the electric field can be written as 



j a (r; u) = J a aP (r, r'; uo) Ep (r' ; uo) dr 1 (62) 
a aP (r, r'; uo) = a D (S a p5{r - r') - V Q V^d (r, r'; u)) (63) 
where d is the diffusion propagator (" diffuson" ) 

V 2 d (r, r'; w) = -5(r ~ r') + ^-d (r, r'; u) (64) 

satisfying the boundary condition 

V„d (r, r'; uo) \ B = V n d (r, r'; u) \ d > = (65) 

Clearly, the linear response described by eqs. ( |62] ) and ( |63| ) is the real-space variant of the non-local current 
response described in Chapter III of Ref. 0). In particular, for an infinite homogeneous system in a longitudinal field 
qE exp (iq ■ r), it is in agreement with the non-local (g-dependent) conductivity given by eq. (3.145) in Ref. Q. It is 
also equivalent to the "generalized Einstein equation" (see eq. ((2^) and below), the latter being its differential form. 

Using eqs. (0), (ph, ( p4| ) and ([35]) and E = — V0, we find from eq. ( |62| ) via integration by parts, 

j (r; u) = -iuoe 2 ^V / d (r, r'; to) <j> (r'; u) dr' (66) 
dfi J 

and 
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p(r;u) = / U(r,r';uj)(j)(r';Lu)dr' (67) 



n(r,r'; W ) = e 2 ^- (-S(v r') + ^-d (r, r'; W )) (68) 

Together with the Poisson's equation, V 2 cf> = — 47rp, eq. d67j ) is an integro-differential equation with respect to cf>. At 
small frequencies, u> -C D / L 2 , it can be solved by successive iterations since, in respective orders of ui, it reduces to 
differential equations (this point is further elaborated in Appendix B since it is related to derivations in the quantum 
limit in Refs. § - §). 

It follows from eqs. (|l3|), @ and (||) that 

D a (r; to) = J e al3 (r, r'; w) E fj (r'; w) dr' (69) 

4:7rz 

e a/3 (r,r';w) = 5 Q/3 5(r - r') a a p (r, r'; u) (70) 



where e~ is the nonlocal dielectric tensor, which satisfies the boundary conditions from (|65|). The electric displacement 
vector D outside the sample coincides with the electric field, D = E°" , and D n is continuous through the boundary. 

For ellipsoidal geometries (see 0), D a =a a pE j3 (as shown above, a a3 — S a s and a a p — 8 a p (l + for a slab and 

a sphere respectively). Then it is possible to invert eq. (|69|) 



E a (r; u) = const (J^tfi ( r > r 'i u ) dr '^j E op 



so that one can find the field inside the sample as a response to external field. A straightforward calculation for a slab 
(see Appendix C) confirms that the nonlocal-? result coincides with that obtained in the phenomenological derivation. 

The correspondence between the integral equations (p3)-(p3) and the phenomenological description in Section II is 
easily shown by applying the V-operator to both sides of eq. ( |67| ) and using eq. (|6^). This gives the "generalized 
Einstein transport equation" ( p3| ) 

j (r; lj) = cr D E (r; w) - D D Vp (r; w) 

while the boundary condition ( |65|) is equivalent to j n \d — 0. Clearly, our response function formalism correctly 
reproduces the combined effect of Maxwell's equation (screening), continuity, and diffusive approximation. 
The diffuson can be written as 

i^i^E^ 1 (71) 

m=0 m D n 

in terms of the volume- normalized eigenfunctions f2 OT) such that 

V 2 n m (r) + X 2 m Q m (r) = (72) 
V n O ro (r) | a = (73) 

We remark in passing on zero mode, Xq — 0, for which f2o = V^ 1 ^ 2 . This introduces a term 



V- 1 J cj)(r-uj)dr 



in the expression of p versus <j), which can be always chosen zero. Dropping the zero mode replaces the diffusion 
propagator by 

, / / \ ( r ) ( r ') 
d(r,r = > ■■ 

The modified diffusion propagator satisfies the following equation 

V 2 d (r, r'; w) = -6(r - r') + i + ^-d (r, r'; lo) (74) 

V L>D 

and is the Fourier transform of the time-dependent diffuson of Ref. || . 

It must be emphasized that the response-function formalism is, obviously, subject to the same condition of appli- 
cability - an assumption of diffusion approximation - as the phenomenological approach. 
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IV. DISCUSSION 



We have obtained closed form expressions for the absorption of the quasi-static longitudinal electric field in a slab 
and sphere. We have also derived corresponding real space expressions for the nonlocal linear response functions. It 
is of interest to compare our present results with the absorption in a transverse field. As pointed out in Ref . [EI , in a 
metal the magnetic component of the transverse field is much larger than the electric field. Therefore, we will address 
the magnetic dipole absorption. 

It can be shown that the magnetic moment M. of a slab and a sphere (!]] in a field H y = Hoe luJt are given respectively 

by 

VH ( tanh (Lx/2) 
M = — 1 - 



8tt V Lx/2 
and 

M = ( 1 _ 3 coth(ax) ( 3 



8tt y ' ax ' (ax) 2 J 
where 

1+i , c 
x = — - — , o = - 

^2-KOoLO 

and 5 is the penetration (skin) depth. Functionally, these results are virtually identical to the ones for the longitudinal 
filed. The main differences are the absence of the depolarization factors and the purely dynamical nature of S. At 
low frequency the latter can be much larger than the sample size, e.g. ax -C 1, in contrast with LA 3> 1 for a good 
metal in a longitudinal field. This means that the transverse field can easily penetrate the metal and, as a result, the 
magnetic dipole absorption 

Qm = ^Rc{iuMH } = -^LuH lm{M} 
can become dominant Indeed, for a spherical metal particle in a wave, Hq — Eq, we find that jjj 

) M tt 2 ( 2a£ 



Qrd 90 VApA" 1 



where A p is the plasma wave-length |1(J. Assuming lot ~ .05, so that 6 — \ v j \J2-k 2 lot ~ \ p where \ p ~ 3000 A 
fiof , and taking A -1 ~ 1 A [0 and £ ~ 100 A (mean-free-path is typically smaller in small particles than in bulk 

o 

materials), we find Qm/Qrd ~ 10 for a ~ 500 A ■ 
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APPENDIX A: ABSORPTION IN TERMS OF ELECTRIC MOMENT 

Since there is no charge outside the sample, the integral in eq. ( |l8|) can be extended to the entire space. We can 
write 

P (t;uj) = -V 2 ${v;lo) 

where the potential 

0(r;w) = 4>(r;u)-(j)o (r;w) 
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is such that (f> (r — > oo; ui) = 0. Using these equations, we find 

Q = Re|-y J p(r; W )^o(r)rfr|+Re|y J V 2 (r; w) 4>* (r; w) drj (Al) 



By the Green's formula, we find 



RG {y/ V2< ^( r ; w )^* ( r ;^) dr | =Re {y / V 2 ?* (r;w)0(r;w)dr| 



= -ReHyy V 2 0(r; W )0*(r; W )drj | = - Re | J V 2 (r; w) 0* (r; w) dr 
which shows that the second term in eq. (|Al|) is zero and proves eq. (|l9|). 

APPENDIX B: SMALL-FREQUENCY EXPANSION 

In the zero's order in u), we find from eq. ) 

Pst (r) = -e 2 ^<P st (r) (Bl) 

and 

V 2 st (r) - A 2 st (r) = (B2) 
This must be supplemented by electric field matching at the boundary. Obviously, j sf = <7o~E st — _DV ' p st = 0. In the 



next order in u, 4> — 4> s t + 4>u> and p — p st + p w , with p w oc w. Using eqs. (p7f) and (B2), we find 



Puj (r; u) = -e 2 ^ [<j> u (r) - cj> RD (r; «)] (B3) 

and the differential equation for W 

V 2 ^ (r) - A 2 w (r) = -K 2 4> RD (r; w) (B4) 

which is easily solvable for simple geometries, such as slab and sphere, and which produces results corresponding to 
the small-frequency expansions of the respective complete solutions. 
Above, we introduced the notation 

fan (r; u) = ^ J d (r, r'; 0) st (r') dr' (B5) 

which we call the Rayleigh-Drude potential. It gives the dynamical field in the bulk and is responsible for absorption. 
As follows from eq. (|65|), this field falls off to zero at the boundary, ^RD,n\d = 0. In a slab, for instance, 

( / )RD (z;u;) = -^((z-^E + ( f> st (z)^ (B6) 

where we used 

oo cos (nip) cos (ni^:) 1 , 2 j s 

« ») - £ 1 1 = " " " + H < B7 > 

m=i l~rJ v 7 



Comparing with expression ( |40D , we observe that eq. ( JBq ) is obtained by expanding the former in oj with the account 
for the w-dependence in ejj, while ignoring such dependence in A (see eq. (p6|)). 
Expanding eq. ( |l8| ) for small frequencies, we observe that 

Q = ^ j-y / [ Pu) (r) st (r) - p st (r) ^ (r)] dr 
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With the use of eqs. (Bl) and (B3), the latter reduces to 
Q = Re | y J p st (r) <j) RD (r; w) drj 



6 sf (r) d (r, r'; 0) <j> st (r') dr'dr =— / / p st (r) d (r, r'; 0) p st (r') dr'dr 



(B8) 

which should be compared with expressions used in Refs. - Qf| (notice also that <fi s t (z) = (—Egz) a (z) for a 
slab and <p s t (r) = (— Eo ■ r) / (r) for a sphere where the functions g and / are easily found from eqs. [ll] and |54| 
respectively). A quick check for a slab confirms that, using eqs. @ and @, eq. @ yields eq. © for M > 1. 



APPENDIX C: DIELECTRIC TENSOR IN A SLAB 



For a slab, 



whereof 



and 



e(z,z';uj) 



e~ 1 (z,z';uj) 



- 2 cos(^)cos(^) 
d(z,z ;wj = > -2 

m=l \ L J ' 1 Do 



E 2 I . n 1 \ fTYVKZ\ I mnz' 

- 1 + A j sin — — sin 



E 2 / . n 1 \ / m7rz \ 

- l + A 2 5 sm — — sir 



The field inside the slab is given by 



= E 1 ( 1 + A 2 r 



m — odd 



( mnz \ 



which, upon summation, yields eq. (38) 



(CI) 



(C2) 



(C3) 



(C4) 



[1] L.D. Landau and E.M. Lifshitz, Electrodynamics Of Continuous Media (Pergamon Press, New York, 1960). 
[2] C.F. Bohren and D.R. Huffman, Absorption And Scattering Of Light By Small Particles (John Wiley & Sons, New York, 
1983). 

[3] A. A. Lushnikov and A.J. Simonov, Phys. Lett. 44A, 45 (1973). 

[4] D. Pines and P. Nozieres, The Theory Of Quantum Liquids, v. I, pp. 192-195 (W.A. Benjamin, New York, 1966). 

[5] R.A. Serota, J. Yu, and Y.H. Kim, Phys. Rev. B42, 9724 (1990). 

[6] L.P Gor'kov and G.M. Eliashberg, Sov. Phys. JETP, 21, 940 (1965). 

[7] A. A. Lushnikov, V.V. Maksimenko, and A. Ya. Simonov, Sov. Phys. Solid State 20, 292 (1978) and V.V. Maksimenko, A. 

Ya. Simonov, and A. A. Lushnikov, Phys. Stat. Sol. (B) 83, 377 (1977). 
[8] K.B. Efetov, Phys. Rev. Lett. 76 1908 (1996). 
[9] Ya.M. Blanter and A.D. Mirlin, |cond-mat/9705250| (1997). 
[10] N.W. Ashcroft and N.D. Mermin, Solid State Physics (Saunders College Publishing, Fort Worth,1976) 
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